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EFFECTS OF A WAVY NEUTRAL SHEET ON COSMIC RAY ANISOTROPIES 

.... J : K ®ta U» 2 ) and J.R. Jokipii <1) 
rl T ? e ^ n ^ ver8 *- t y of Arizona, Tucson, Arizona, USA 
(2) Central Research Institute for Physics, Budapest, Hungary 

Afe£tCAS^.We present the first results of our 3-D numerical code calcula- 
inffh ?• ani8 ? tr °P ie *‘ The code includes diffusion, convection, 

C00lln8, and drift in an interplanetary magnetic field model 

all the Sin D ® Utral ? heet * We find that the 3-D model can reproduce 

all the principal observations for a resonable set of parameters. 

d Sf?f U K tiOP ‘ In the la8t . decade, the effects of curvature and gradient 
drifts became a central issue in the theory of cosmic ray transport. It 

n^nt eaa8u ? 8e8ted . that drift ma y Play an important, and perhaps domi- 

resMnsible f J°?h 1C ^ P r °P a 8 ation \» the heliosphere, and it may be 
responsibie for the asymmetries appearing in consecutive 11-year cycles 

Snd Hubbard , 197 7; Jokipii and Kopriva 1980; Kota 1979; 

Lift" model! X?", 19815 K V &ad J ° ki P Li 1983 >‘ Tbe first success of 
tion of i-h t P ex Pl ainin 8 galactic cosmic-ray phenomena was the explana- 
• . f the Phase shift of the solar daily variation (Levy 1976) observed 

in the years of the seventies, following the polarity reversal of the 
solar magnetic field. The first quantitative 2-dimensional (2-D) calcula- 
tion was carried out by Jokipii and Kopriva (1980). In this work, howe- 
ver, the too small value of the diffusion coefficient, rc, led to unrea- 

fo^e-fi* n i l / 0 »n r i 0P /* 8 in . same ca8es * Kota < 19 81) derived an approximate 
1- ! ° n v Wlth 8 virtuall y perfect isotropy. This model, 
however, relied upon the too simple picture of 'hard-sphere' scattering. 

Pomeranfc 611 Jf Phase shift of the solar daily variation (Duggal and 
Pomerantz, 1975) is naturally explained by 2-D models (Levy 1976, 

u\o1 19U l magnetic -“fixation of the seventies 
eonvi^i i d K \r all r radial density gradient which cannot balance the 
AnTw i? -lar wind and thus results in a net outward streaming. 

other well established observation is the presence of the polarity 
dependent N-S anisotropy associated with the B x Jn streaming 

Sirr tCh 1970 ’ Pomerantz and Bieber 1984). In the seventies (A>0), 
thi! 8tr . eaming “ directed away from the neutral sheet. In a 2-D model, 
this pattern of streaming is hard to reconcile with the div S < 0 

requirement, thus 2-D models are bound to encounter difficulties in 
explaining both observations. 11 in 

ting h th! a «ti«! iff i C w 1 c y * in Pri^ipfe at least, may be removed if, viola- 
ting the axial and N-S symmetries, a wavy neutral sheet is included. It 

aM "to PUrP °r ° f J? 18 W ° rk tC demon8trate that a 3-D model is indeed 
able to reproduce all components of the observed anisotropies. We present 

sheet I' 8 1 ani80tr °Py re8ult8 °f °ur 3-D code incorporating" a Cavy neutral 

U8ad * 3 " D nu “ erical cod ® to 8 ®lve the modulation equation 
including diffusion, convection, adiabatic deceleration and drift. The 

?°£ al and T tke . 8C . h ®“ e 0 , of calculation were described in detail elsewhere 

issuel Brieflt" 11 1983 ? *** a , l8 5’. the Preceding paper SH-4.2-10 in this 
a «i y * 8 “ 8Ual 8 P xra ,l ^ le l^ ia adopted, the magnetic equator is 

J£ ed pla “ at the 8un » which then evolves into a wavy sheet (Jokipii 
Thomas 1981). The case of A>0 corresponds to outward polarity above 
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the sheet and inward polarity below the sheet while A<0 corresponds to 
the opposite configuration (sixties and eighties). Steady state is 
assumed in the frame corotating with the sun. The most serious limitation 
of the code is that it assumes constant solar wind speed thus many 
phenomena, like shocks, are precluded. 

Calculations were carried out for protons in the 1 - 10 GV range. The 
parallel diffusion coefficient, tc H , was assumed to be inversely 
proportional to the magnetic field strength, B, 

- Ko pl/2 P <B e arth/ B > 

with P being the particle rigidity in GV, p is the particle velocity 
in units of velocity of lights and K 0 is a normalization constant in 
the range of lO 2 * - 10 2 ^ cm 2 /sec. The ratio of the perpendicular and 
parallel diffusion coefficients was kept constant at “0.05-0.20. 


Result 8 and Discussion. Anisotropies were calculated at three 
heliocentric distances (0.5, 1 , and 5 AU) over the full range of 

heliographic latitudes and longitudes. Here, we present the results near 
the earth. The anisotropies to be reported are obtained at the 
helioequator, at 1 AU, and averaged over longitudes in a magnetic sector. 

The ecliptic components of the anisotropy responsible for the solar 
daily variation are given in Figures 1 and 2 (P*=2.3 GV; “.05). It 

should be noted that the anisotropies obtained for a flat sheet (dashed 
lines) show sharp changes at the neutral sheet. The actual values (dots) 
may considerably differ from the averages over a (-5°, 5°) latitude band 
(open circles). In most cases, A>0 gives an earlier phase and a smaller 
amplitude which is in general agreement with observations. Similar re- 
sults were obtained for other rigidities, too. At large values of K Q> 
understandably, drift effects diminish and a near perfect corotation 
applies. The breaks in the lines in the K q = 1.5-5.10 22 cm 2 /sec range 
indicate that corotafion should be reached somewhere in this interval. 

Figure 3 shows , the zenith angle component of the average near 
earth anisotropy above the neutral sheet for P “ 2.3 GV. Consider first 
the case A>0, vhen the ’observed 1 value of is negative in accordance 
with the sense of the B. xVn streaming. Curve^(a) corresponding to 
« .05 and ot * 15° yields the correct sign for£<*. Larger tilt angle ( oc » 
30°, curve (b)), however, may already give positive values, too. If we 
take, on the other hand, a larger perpendicular diffusion (curve (c): 
K, x /«Si 20, ot =30° ),£<$, will agaih point in the proper direction for all 
values of K 0 . The underlying physical picture is that, in the case of 
large tilt angles and small perpendicular diffusion, most particles reach 
the earth without having interacted with the neutral sheet. Being too 
far from the earth, the neutral sheet becomes irrelevant for A>0. A6 for 
A<0, particles intersect the sheet, several times before reaching the 
earth. As a result, the calculated always shows the proper sense and 
is fairly independent of the tilt angle, ©c . 

Figure 4 indicates that the magnitude of£^ increases and its sign 
becomes more distinctive with increasing rigidities. This finding can 
also be anticipated since drift effects are expected to be more 


pronounced at higher rigidities. 

The typical azimuthal dependence of the N-S anisotropy is presented 
in Figure 5. In most cases, we cannot find a one-to-one correspondence 
between the sign of and the polarity of the field. In general, smal- 
ler tilt angle, and larger results in a better correlation; the 
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correlation also improves at higher rigidities. At 5 AD heliocentric 

“‘T:\J irtUally 811 ° Ur runs gave 8 100 P er cent correlation. To 

Svel7 e tiahter ?t ¥** tha \. the waves in the neutral sheet become rela- 
^J*V at larger distances from the sun. Thus we expect the 

effects of waviness to be more direct there. 

C onclusion. Our numerical results demonstrate that the inclusion of a 

Dies T?- 1 j 11661 may ex P lain «H components of the observed anisotro- 
I^WrA*/ 8eneral agreement between 'theoretical* and 'observed' 
above^K i f " a V de range of Parameters. At high latitudes, well 
toward H? el0W he 8l ^ eet » how ever, we predict the anisotropy to point 
n , aid / e equator in the case of A>0. This is in contrast to the 
poleward direction expected from the B x Jn term only. 
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Figure 1. Radial and azimuthal components of the average near earth 

“J!rr P fi C ! 1C , Ulated f ? r A>0 ’ P “ 2 * 3 GV » 'V^-OS. Dashed and solid lines 
f®2? t 9, flat (see text) and way y sheets, respectively. K i 8 in unitfi of 
10 cm /sec. The phases of the resulting daily waves are also shown. 
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Figure 2. Same as Fig. l.» for A<0. 
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Figure 3. The average value 
of 1^. for the earth being 
above the neutral sheet. 

P-2.3 GV, K e « i“ unit8 
of 10 22 cm 2 / sec. 
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Figure 4. vs rigidity cal- 
culated for K o -10 22 cm 2 / sec, 

Kj. /ic*«0.05. 


Figure 5. Azimuthal dependence of 
calculated for P-2.3 GV, ic^/iCj-.OS. 
(a) AX), o6-l 5° ; Cb) AX>,oc«30°; 
(c) A<0,ot -15°. 










